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文章 编号 :1005-3085(2010)05-0820-07 
埃 尔 米 特 广义 汉密尔顿 矩阵 的 广义 逆 特 征 值 问题 * 


魏 平 !， 张 忠志 2， 谢 冬 秀 3 
(1- 华南 理工 大 学 数学 系 ， 广 州 510640; 2- 东莞 理工 学 院 计算 机 学 院 ， 东 莞 523808; 
3- 北京 信息 科技 大 学 理学 院 ， 北 京 100192) 
摘 要 : 本 文 利用 埃 尔 米 特 广义 汉密尔顿 矩阵 的 性 质 与 矩阵 的 分 解 理论 ， 导 出 了 埃 尔 米 特 广义 汉密尔顿 矩 
阵 的 广义 道 特征 值 问题 解 的 一 般 表 达 式 。 进 而 运用 希 尔 伯 特 空间 的 通 近 理论 ， 对 任意 给 定 的 于 阶 
复 和 矩阵 对 ， 证 明 相关 最 佳 逼 近 解 的 存在 性 与 惟一 性 ， 得 到 了 最 佳 逼近 解 的 表达 式 。 
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矩阵 的 道 特征 值 问题 与 广义 逆 特 征 值 问题 在 工程 技术 中 有 着 广泛 的 应 用 ， 关 于 这 类 问题 的 
研究 已 取得 一 系列 成 果 。 例 如 文献 [1,2] 分 别 就 双 对 称 和 矩阵 和 对 称 正 交 对 称 和 矩阵 的 逆 特 征 值 问题 
进行 了 研究 ， 给 出 了 可 解 的 充分 必要 条 件 。 文 献 [3,4] 利用 埃 尔 米 特 广义 汉密尔顿 矩阵 的 结构 ， 
分 别 就 埃 尔 米 特 广义 汉密尔顿 矩阵 逆 特 征 值 问题 的 可 解 条 件 和 线性 矩阵 方程 的 埃 尔 米 特 广义 汉 
密 尔 顿 矩阵 解 进行 了 研究 ， 并 得 到 了 它们 的 最 佳 逼近 解 。 而 文献 [5,6] 则 分 别 就 实 对 称 和 矩阵 的 广 
义 逆 特 征 值 问题 和 反对 称 正 交 反对 称 和 矩阵 的 广义 逆 特 征 值 问题 进行 了 研究 ， 推 得 了 问题 解 的 一 
般 表达 式 和 最 佳 逼 近 解 。 本 文 就 埃 尔 米 特 广义 汉密尔顿 矩阵 的 广义 逆 特 征 值 问题 进行 探讨 。 

H Rg"*", Cnxm 分 别 表示 所 有 m x mr k EBE. SÓBPERISEG. ERFARA E 
和 Moore-Penrose 广 义 逆 分 别 记 为 4 和 A+， 表 示 n 阶 单位 矩阵 。UC"x"* 表 示 所 有 nn 阶 
西 和 矩阵 的 集合 ，rank(4) 表示 和 矩阵 4 的 秩 ， 对 任意 矩阵 4 c Or, B c C, AS B 
zz A Ej B ËJ Kronecker 积 


AQB = (aB) e C"?*", i=l, m, j=1,.,n, V A,B 8 Onxm, 


(A, B) = tr(BP A) 表示 4 与 BB 的 内 积 ， 不 难 证 明 由 此 内 积 诱导 的 矩阵 范 数 4|| = /(A, A) 是 
Frobenius 范 数 ， 且 在 该 范 数 下 C"*x™ 是 一 个 希 尔 伯 特 空间 。 令 OASR"*" 表示 n 阶 实 正 交 反 
对 称 和 矩阵 集合 ， 即 


OASR""" = {J| JTJ = JJT = In, J = —J77, J e R^*"). 


定义 1 给 定 Je O4SRnxn， 若 4E Onn il AE = A, JAJ = A#， 则 称 A 为 n 阶 埃 尔 
KE UE RERE, BD nEBAZORNH XDURARSDEEERIS MU HHO”, 
定义 2 WA- (045), 7 ido; = (a1, 025 ° ,Qni), i= 1,::- ,m, 令 


vec(A) = (o1, 02,77, am)”, 
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则 称 vec(A) 为 矩阵 4 的 列 拉 直 。 

定义 3 给 定 JeO45SRn"xnmn， 对 任意 的 ze Cn, i 

1) #iJz= xz， 则 称 z 为 J- 对 称 向 量 ，n 维 对 称 向 量 的 全 体 记 为 SC”"。 

2) 车 iJz = 一 xz， 则 称 z 为 J- 反 对 称 向 量 ，n 维 反对 称 向 量 的 全 体 记 为 ASC" , 

TER, RS HHO” EO 的 一 个 线性 子 空 间 ， 而 且 它 与 矩阵 J 有 关 ， 全 文 总 是 假定 矩 
阵 了 是 不 变 的 。 运 用 矩阵 了 的 性 质 ， 有 .有 = —L,, Wik, nd BER, 

4 

P, = i), Pa = i-i) 


其 中 i = /-1, 

不 难 证 明 号 和 已 都 是 正 交 投 影 矩 阵 ， MEA P, + P; = I, P, P, = 0。 因 此 ， 存 在 单位 
FERIERE UL, Uo € Cnxk (n = 2k), 使 得 忆 = U,UP, P, = U3UT , WU = (U1, U2). WE 
ÉE Pi, 忆 的 性 质 ， 容 易 证 明 U 是 nn 阶 西 和 矩阵 ， 即 UE UC"*", 

问题 1 给 定 


X = (z1, £2, sata Cm) € Cm. 人 一 diag( 和 il, As, Am) € Qmm 


XABPE A, Be HHO?*" 使 得 
AX = BXA. (1) 


问题 2 ”对 于 任意 给 定 的 矩阵 A, Becr, cR (A, B) € SA p 使 得 


(A.B) - (4, B) = min — (4 B) - 4, B)]| (2) 


其 中 54.B 是 问题 1 的 解 。 

本 文 结构 如 下 : 在 第 二 节 利 用 埃 尔 米 特 广义 汉密尔顿 矩阵 的 性 质 推 导 问 题 1 解 的 一 般 表 达 
式 。 第 三 节 证 明 问 题 2 解 的 存在 性 与 惟一 性 ， 并 给 出 最 佳 逼 近 解 的 表达 式 。 第 四 节 给 出 求解 问 
题 2 的 算法 及 数值 例子 。 

2 问题 1 的 解 


引 理 1 设 Ae Omxn be€ Cm， 则 线性 方程 组 Az = 有 解 的 充分 必要 条 件 是 44+0 = 
b。 这 时 方程 组 的 通 解 为 z = Atb+ (I — 4+4)z， 对 任意 的 ze C". 


引 理 2l8] 24374 
0 N 
S-(A|A-U UE NUES. (3) 
NH 0 
WAS = HHOnxn, 


不 难 证 明 下 面 引 理 。 
引 理 3 已 知 W eC"*", Rec?*, WIT- RJ? - ||T — W|? = min, T e C"* 的 解 为 


T- 3(R+ W). (4) 


33:84 ER A, Be HHOnxn, WIERE AURE AX = BXA， 其 特征 向 量 可 
表示 为 J- 对 称 向 量 或 J- 反 对 称 向 量 。 
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证 明 H Az = ABr Kk A, B c HHO?*" 41 JAJz = XJBJz， 进 而 4(Jz) = AB(Jz), HJ 
此 ， 车 z 是 4, B 的 广义 特征 向 量 ， 则 恒 有 .Jz 也 是 4, B 的 广义 特征 向 量 ， 所 以 zx 与 Jz 必定 线 
性 相关 。 则 存在 不 全 为 零 的 复数 k1, ko 使 得 kiz 十 k2Jzx = 0， 两 边 同 时 左 乘 以 J] 得 ki1Jx 一 kzx = 
0， 由 这 两 式 可 得 如 十 上 把 = 0， 因 此 ，k1 = iko WA iJe = +z, 


2 


MILI 
一 以 > 人 


X jf o) QF = HEQE, (5) 
—X3^ 0 0 


H = (Hı, H, ) € Ug. Q = (Qi, Q2 ) € UOTA HË = (S1, S2) € Qo-n)xn 
Tj; Tm—T|i 


Ti T7L—Ti 


其 中 


再 设 
r2 = rank((ST & XP) — SZ & (AF xXP)), 


(ST & XF) — ST e (AF XE) 的 奇异 值 分 解 为 


T 0 
(ST & XF) — ST &(AEXP) =E Ë ') FP = ETFH, 


其 中 
F- (F, F> ) € Qr n-ni)xk(n-n). 
T2 k(n—ri)—r2 
定理 1 设 
X —(z1,22,5::, 24) E C”*™, A = diag(à1, MN , Àm) € Cn, 
则 问题 1 的 解 为 


"En P Pus Sp m. ^ “es vn. (6) 
(GS1)” 0 (G52) 0 


其 中 Ge ckx) 满足 vec(G) = FK, K € Ceon- 为 任意 复 列 向 量 。 


第 5 期 魏 平 等 ， 埃 尔 米 特 广义 汉 密 尔 顿 矩 阵 的 广义 逆 特 征 值 问题 823 





证 明 由 引 理 2 知 ，4, B 可 表示 为 


H AX = BXA 有 


MJ 


即 有 
N. X> = N3X2A, (8) 
NP X1 = NË XA. (9) 
H (8) 式 得 
X 
(N1, No) ( i ) =0, 
—X;A 
所 以 
(Ni, N.) = GH = G(S,, S2), (10) 


其 中 51, S, e C(in-ri)xk G e Ckx(n—ri) 为 任意 矩阵 。 
由 (9) 式 有 XE Ni = AFP XP NS, EK 0) 式 代 入 得 


XHGS! = AF XE GS, (11) 
进而 
[(ST & XP) - ST ® (A? XP)]vec(G) = 0, (12) 
则 由 引 理 1 知 (12) 式 的 解 为 
vec(G) = FaK, KEeOrn-r)-r (13) 


将 (10) 式 与 (13) ARA (7) 式 即 得 问题 1 的 解 (6) Ro 


3 问题 2 的 解 
定理 2 任意 给 定 和 矩阵 A. B e Cr"x"， 则 问题 2 有 唯一 的 最 佳 逼近 解 ， 且 这 个 解 可 表示 为 


A 0 Gs z 0 âs 
A=U| .py “Jun B=U| _ a V QUE (14) 
(Gs) O (GS) O 
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其 中 Ge ch 满足 


F FH [vec((UB AUS, UE BU2) H2) + vec((UP AP US, UB BF U3) H2)] 


vec(G) — > 





(15) 


证 明 ”由 定理 1 和 下 - 范 数 的 正 交 不 变性 有 


| (4， B) m (A, B) 


l4 - 4| + ||B - B| 
0 GS; UHAU! UFAU, " 0 GS; UBBU, U¥ BU, 
(Gs 0 U# ĀU, U# ĀU (GS) 0 U# BU, U# BU 


||GS1 - u# Aus| + ||GS2 — UF Bus + GS — UF A" U|? + ||Gs; — u# B” Ul 

















HUF Au |" + |o Ava? + uf Bon]? + jv? Bual 

lIG(S1, 82) — (UË Aus, UF BU2)||? + ||G(S1, S2) — (UF A” Ua, UF BF Uz)||" 

«Ju Av. + os Āu? + U Bus f + [|y Bul 

|GHZ H — (UF AUS, UF BU H ||? + ||GHZ! H — (UË AF Us, UF B"Us)H|| 
«Ju? Av. |f + j£ Ava]? + UF Bus [f + jo Bual? 

||G — (UF Ava, UF Bv) n; | + || — (Uf AP us, UF B" uz) Hal)? 

十 | (UF iuz, UF Bus) m | + (uf AP us, UF B" us) nn ||" 

«Ju? Av? + os A| + jv? Bus | + os Bua 

EK — vec((UF AUs, UF BU2)H2)||? + ||F2K — vec((UF A"Us, UF B” U2) Ha) ||? 
+||(UE Aus, UE Bus) n ||? + |(uF AP va, UE B” uz) |" 

«uf Au, * + of Avo + uf Bo]? ju? Bua 

|K — FF vec((UP AUs, UF BU2) Hz)|| + ||K — Fë vec((UF A" Us, UF B"U;) Ma) ||? 
4| FE vec((UF AUs, UE BUo)Ho)|| || FF vec((Uf A" U2, UF B"Uz)Hz)||" 
+||(U# Ava, v? Bus) | + | (Uf AF us, UF B" U) | 


*|u£ Aui |? + už Au]? + | BU, ||? + ||uz Bull’. 


由 引 理 3 知 PN 
amm, ,|(4, B) - (A B)|, 
ZERA 
g — Fflvec((Uf! AU2, UF BU) Ha) + vec((Uf' A" Us, Uf! BH U2) H2)) (16) 


2 


将 (16) ARA (7) 式 即 得 问题 2 的 解 。 
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4 ”算法 及 数值 例子 


算法 求解 问题 2 的 步骤 如 下 : 
输入 J, A, B, X, A; 
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2 HP, = )MI+iJ) P, = 3M(I-—iJ), W P, P>; 
3) P, =U,UE, P, = UsUE, U = (U1, U2) VPE U,, U2, U; 
4) 由 (5) 式 计算 H5, $1, S9; 
5) 由 (15) 式 计算 G; 
6) 最 后 按 (14) 式 计算 A, B. 
例子 n=4, m= 2, 上 二 2， 给 定 
0 0 1 0.7071 — 0.7071 
0 10 —0.7071  —0.7071 2+i 0 
人 一 , X= s , 
—1 0 0 0.7071i —0.7071i 0 3-8i 
-1 0 00 —0.7071i 0.7071i 
0.5310 + 0.4394i 0.3697 + 0.1192i 0.8380 + 0.4301i 0.4213 十 0.8485i 
i 0.7628 + 0.1702i 0.7785 + 0.6351i 0.5283 + 0.99761 0.6503 + 0.4354i 
0.9137 + 0.1837i 0.9785 + 0.07701 0.2035 + 0.3861i 0.9605 + 0.47401 | ` 
0.7174 + 0.0825i 0.0578 + 0.1941i 0.0660 + 0.5052i 0.3904 + 0.5693i 
0.1572 + 0.7103i 0.5637 + 0.0807i 0.0781 + 0.5058i 0.0432 + 0.7366i 
5 0.0791 + 0.9287i 0.6270 + 0.2687i 0.5883 + 0.6907i 0.5444 + 0.7839i 
0.0457 + 0.1103i 0.3349 + 0.2047i 0.7257 + 0.4955i 0.6301 + 0.9093i 
0.2786 十 0.0951i 0.4144 + 0.9163i 0.2733 + 0.6343i 0.1257 + 0.3191i 
按 以 上 算法 得 结果 如 下 
一 0.2887 0.9353 0.1293 — 0.0647i 0.1293 — 0.0647i 
Sı = | 0.5774 + 0.2887i 0.1293 + 0.06471 |, S> = 0.6767 —0.3233 
0.5774 + 0.2887i 0.1293 + 0.0647i —0.3233 0.6767 
—0.1841 — 0.4505i 
G= 


0.2543 — 0.9033i 


0.1647 十 0.5271i 


0.0772 — 0.0693i 


D + 0.7948i 


最 后 得 A, B 的 唯一 最 佳 逼近 解 A, B 28 


—0.0341 0.1027 — 0.0433i 
ros = + 0.0433i 0.2395 
0.6382 — 0.0368i 0.7700 
0.5063 0.6382 + 0.0368i 


—0.1027 + 0.0433i 


—0.1077 + m 


0.6382 4- 0.0368i 0.5063 
0.7700 0.6382 — 0.03681 
—0.2395 —0.1027 — 0.0433i 


0.0341 
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一 0.0554 一 0.0410 — 0.0451i 0.2909 十 0.0200i 0.2313 
s= —0.0410 + 0.0451i —0.0266 0.3506 0.2909 — 0.0200i 
0.2909 — 0.0200i 0.3506 0.0266 0.0410 — 0.0451i 
0.2313 0.2909 + 0.0200i — 0.0410 + 0.0451i 0.0554 
参考 文献 : 


[1] HBK, KE WAE. 双 对 称 矩 阵 逆 特征 值 问 题解 存在 的 条 件 [可 ., 计算 数学 ，1998, 20(4): 409-418 
Hu X Y, Zhang L, Xie D X. The solvability conditions for the inverse eigenvalue problem of bisymmetric 
matrices[J]. Mathematica Numerica Sinica, 1998, 20(4): 409-418 

[2] 819936. Ka AER. IRERE AREE RI]. 计算 数学 ，2003, 25(1): 13-22 
Hu X Y, Zhang L, Zhou F Z. The inverse eigenvalue problem of symmetric ortho-symmetric matrices[J]. 
Mathematica Numerica Sinica, 2003, 25(1): 13-22 

[3] Zhang Z Z, Hu X Y, Zhang L. The solvability conditions for the inverse eigenvalue problem of Hermitian- 
generalized Hamiltonian matrices[J]. Inverse Problems, 2002, 18: 1369-1376 

[4] Zhang Z Z, Hu X Y, Zhang L. On the Hermitian-generalized Hamiltonian solutions of linear matrix equa- 
tion[J]. SIAM J Matrix Anal Applm, 2005, 27(1): 294-303 

[5] WE. 谱 约 束 下 实 对 称 矩 阵 东 的 最 佳 逼近 [J]. 高 等 学 校 计算 数学 学 报 ，1990, 2(6): 177-187 
Dai H. Optimal approximation of real symmetric matrix pencil under spectral restriction[J]. Numerical 
Mathematics: a Journal of Chinese Universities, 1990, 2(6): 177-187 

[6] 李 伯 忍 ， 胡 锡 炎 ， 刘 学 杰 . WAR F BOSEBRIE 26 BOSUBIÓB PEOR bJ t Hol [J]. 数值 计算 与 计算 机 应 用 ，2007， 
28(4): 282-289 
Li B R, Hu X Y, Liu X J. Optimal approximation of anti-symmetric and orth-anti-symmetric matrix under 
spectral restriction[J]. Journal on Numerical Mathods and Computer Applications, 2007, 28(4): 282-289 

[7] RE. 矩阵 论 [M]. 北京 : 科学 出 版 社 ，2001 
Dai H. Matrix 'Theory[M]. Beijing: Science Press, 2001 

[8] Zhang Z Z, Hu X Y, Zhang L. Least-squares solutions of inverse problem for Hermitian generalized Hamil- 
tonian matrices[J]. Applied Mathematics Letters, 2004, 17: 303-308 


Generalized Inverse Eigenvalue Problem for Hermitian Generalized 
Hamiltonian Matrices 
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Abstract: By means of properties of Hermitian generalized Hamiltonian matrices and the matrix de- 
composition theory, we derive expressions of the generalized inverse eigenvalue problem for Hermitian 
generalized Hamiltonian matrices. By using the Hilbert space approximation theory, for any given 
square complex matrices, we prove that there is only one optimal approximation solution, and further- 
more derive expressions of the optimal approximation solution. 
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